Band offset of GaAs/Al^Gai-^As heterojunctions from atomistic first principles 
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Using an atomistic first principles approach, we investigate the band offset of the GaAs/Al^Gai-^As het- 
erojunctions for the entire range of the Al doping concentration < x < 1. We apply the coherent potential 
approach to handle the configuration average of Al doping and a recently proposed semi-local exchange poten- 
tial to accurately determine the band gaps of the materials. The calculated band structures of the GaAs, AlAs 
crystals and band gaps of the Al^Gai-^As alloys, are in very good agreement with the experimental results. 
We predict that valence band offset of the GaAs/AUGai-! As heterojunction scales with the Al concentration x 
in a linear fashion as VBO(x) ~ 0.587a;, and the conduction band offset scales with a; in a nonlinear fashion. 
Quantitative comparisons to the corresponding experimental data are made. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.55.Eq, 71.23.An 
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Introduction. The properties of III-V compound semicon- 
ductors and their heterojunctions have been relentlessly inves- 
tigated for several decades due to their wide-ranging applica- 
tions in electronic and optoelectronic technologies. One of 
most important electronic property of heterojunctions is the 
band offset which describes the relative alignment of the elec- 
tronic bands across the junction interfaced An accurate deter- 
mination of the band offset is critical for understanding quan- 
tum transport properties of the heterojuncton. For many III- 
V materials systems, the band offset has been carefully mea- 
sured experimentally i 2 ^ 

Theoretical calculations of band offset have always been a 
serious challenge. This is because first principles method of 
density function theory (DFT) with local-density approxima- 
tion (LDA) 4 and generalized gradient approximation (GGA) 5 
underestimates the band gap (E g ) of semiconductors. With- 
out a correct calculation of E g for individual semiconductors, 
the calculated band offset between two semiconductors may 
be compromised. Advanced methods such as GW 6 and/or hy- 
brid functional^ can yield accurate E g for many systems but 
require significantly more expensive computation. Another 
difficulty is when there are impurities: the predicted physi- 
cal results must be averaged over the multitudes of impurity 
configurations which is extremely costly in computation. Re- 
alistic semiconductors all have impurities and if the impurity 
concentration is small, one needs to compute systems of large 
number of host atoms in order to accommodate a few impu- 
rity atoms. Because of these issues, predicting band offset has 
persisted to be a challenging theoretical problem. 

Considerable theoretical efforts have been devoted in the 
literature to correctly predict E g . Apart from the GW— and 
hybrid functionals 7 , for pure semiconductors the recently pro- 
posed modified Becke-Johnson (MBJ) semilocal exchange 
potential was shown to give quite accurate E g values for many 
compounds with a computational cost similar to that of LDA. 8 
Even though MBJ is not a fundamental solution to the issue 
of electron correlation, it is practically useful for calculating 
Eg. To deal with the prohibitively large computation required 
for performing configuration average of doped semiconduc- 



tors, one wishes to obtain the averaged physical quantity with- 
out individually computing each impurity configuration by 
brute force as in the super-cell approach. In this regard, a 
widely used technique is the coherent potential approximation 
(CPA)2 as implemented in Korringa-Kohn-RostokeriS or lin- 
ear muffin-tin orbital (LMTO>i± DFT methods. CPA is a sta- 
tistical effective medium approach such that an atomic site has 
x% chance to be an impurity (e.g. dopant) atom and (1 — x)% 
chance to be a host atom, the configuration average is carried 
out analytically hence disorder effect can be calculated for any 
concentration x. Very recently, Ref. |T2| has combined CPA 
with MBJ and reported the calculation of E g for the semicon- 
ductor In-cGai-xN, the results are in excellent agreement with 
the measured data for the entire range of < x < 1. 

In this paper, we employ the CPA-MBJ approach to quanti- 
tatively calculate the band offsets of two semiconductors with 
impurity doping. In particular, we consider the most impor- 
tant heterojucntion, between GaAs and ALGai-^As, and cal- 
culate the band offset as a function of the Al concentration x. 
Even though this heterojunction has been the subject of ex- 
tensive past investigations, theoretically the calculation was 
usually done at a few special values of x where the super-cell 
approach could be applied. Here, the CPA-MBJ approach al- 
lows us to determine a continuous curve of the band offset 
for the entire range of < x < 1, which has not been de- 
rived before from atomic first principles. Our calculated E g 
of ALGai-zAs and the calculated band offsets of the hetero- 
junctions for the entire x range, are quantitatively and excel- 
lently compared with the corresponding experimental data. 

Calculation details. The lattice constants of the zinc-blende 
AlAs and GaAs were individually relaxed by a projector aug- 
mented wave method with LDA as implemented in the VASP 
software^ and the optimized lattice constant of 5.63 A for 
GaAs were used for all the compounds in our calculations 
since the optimized lattice constant of AlAs (5.64 A) is very 
close to that of GaAs. To calculate the band offset with Al 
doping, we apply the CPA-MBJ approach as implemented in 
the Nanodsim software package-^ where the DFT is carried 
out within the TB-LMTO scheme under the atomic sphere ap- 
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these materials. It should be noted from Fig.l(a,b) that the 
band gaps were underestimated by LDA. Next, we apply the 
MBJ functional 8 *^ to calculate the electronic structure agaii>2£ 
and the MBJ restuls are shown in Fig. l(c,d). Compare with 
the LDA bands in Fig.l(a,b), the opening of band gap is ev- 
ident. From Fig.l(c,d), the band gap values at T, X and L 
points are listed in Table Q] The MBJ band gaps are in good 
agreement with the experimental values at the J? and X points, 
and within 7% for GaAs and 15% for AlAs to the experimen- 
tal values at the L point. In contrast, the LDA gaps typically 
underestimate the values quite significantly which is a well 
known issue of LDA. 



FIG. 1: (color online) (a,b) The band structures obtained with LDA: 
(a) for GaAs, (b) for AlAs. Red line is obtained by VASP, blue dots 
obtained by Nanodsim. (c,d) The band structures calculated with 
MBJ by Nanodsim: (c) for GaAs and (d) for AlAs. 



proximation (ASA) Ji For technical details of the Nanodsim 
algorithm we refer interested readers to Ref. Il4lll5l 

In the DFT-CPA-MBJ calculations^, the primitive cell of 
the zinc-blende structure was used to calculate the band struc- 
tures and Eg. To determine the band offset, a (110) system 
containing 9 layers of GaAs and 9 layers of ALGai-^As 
are used to calculate the potential profile through the 
heterojunctionii A12x 12x12 fc-mesh and a 12 x 12 x 1 
fc-mesh were used to sample the Brillouin Zone for the prim- 
itive cell (bulk) and the heterojunction, respectively. For the 
ASA, vacancy spheres were placed at appropriate locations 12 
for space filling, and the same radius were used for all the va- 
cancy spheres and atomic spheres. Spin orbital coupling was 
not considered in this work. 

A technical issue of MBJ semilocal exchange potential is 
worth mentioning. MBJ potential has the following form 8 , 



.MBJ i 



V- J (r)=cv^(r) + (3c-2) 




(1) 



where subscript a is spin index, p a is the electron density, t a 
is the kinetic energy density, and v^(r) is the Becke-Roussel 
potential. 1 - The relative weight of the two terms is given by 
a parameter c which depends linearly on the square root of 
\Vp\/p. For all solids investigated by MBJ potential so far^ 
it appears that E g increases monotonically with c. The c pa- 
rameter can be determined self-consistently as discussed in 
Refs. 1811191 but in this work we fixed its value to reproduce 
the experimental band gaps of GaAs, AlAs and, afterward, we 
used the same value to calculate the electronic structures of 
the alloy ALGai^As and the band offset of the heterojunc- 
tions. 

Band structures of GaAs and AlAs. We begin by calcu- 
lating the band structures of GaAs and AlAs with LDA using 
both VASP 13 and Nanodsim 14 electronic packages. Fig.l (a,b) 
show a perfect agreement of the valence bands and a good 
match of the conduction bands between these methods, con- 
firming that the ASA employed in the TB-LMTO approach 
(Nanodsim) is accurate for calculating physical properties of 



TABLE I: Energies of the conduction band minima at the F, X, and 
L points with respect to the valence band maximum at the T point in 
units of eV, calculated by DFT with the LDA and MBJ functionals 
at zero temperature. The column of LDA V were obtained by VASP, 
other results were by Nanodsim within the TB-LMTO implementa- 
tion. The last column are the experimental values.— 



Material 


E 3 


LDA" 


LDA 


MBJ 


Expt. 2 


GaAs 


r 

X 

L 


0.493 
1.334 
0.948 


0.761 
1.346 
1.100 


1.518 
1.960 
1.691 


1.519 
1.981 
1.815 


AlAs 


r 

X 
L 


2.014 
1.312 
2.086 


2.300 
1.307 
2.191 


3.099 

2.258 
2.835 


3.099 
2.24 
2.46 



Band gaps of Al x Ga\- x As. Having accurately determined 
the band structures for the pure materials GaAs and AlAs with 
the MBJ functional, we apply the CPA-MBJ approach dis- 
cussed above to calculate the the electronic properties of the 
alloy ALGai-^As. The calculated band gaps are plotted ver- 
sus the Al concentration x in Fig. 2. It is important to recall 
- as shown in Fig.l, that GaAs is a direct-gap semiconductor 
where the T valley is lower than the X valley in the conduc- 
tion band; AlAs is an indirect-gap semiconductor where the 
opposite is true, namely X valley is lower than the V valley. 
By increasing the Al concentration x, the alloy changes from 
direct-gap to indirect-gap at a crossover point where the con- 
duction band minima at T and X points have the same energy. 
The existence of the crossover point, at near x ra 0.36, is 
quite evident in Fig. 2. Such a crossover of band gap behav- 
ior is actually known experimentally and, as stated in Ref. y, 
the experimental data implied the T-X crossover composition 
to be at a; = 0.38 at low temperature and at x = 0.39 at 
300 K. It was also experimentally known that the E g scales 
linearly with x in the direct-gap regime and quadratically in 
the indirect-gap regime^ 2 -!- Fitting the calculated results of 
x < 0.3 by a linear function and x > 0.6 by a quadratic func- 
tion, as shown by the solid lines in Fig. 2 excellent consistency 
to these scalings is obtained. In particular, the data points 
of x = 0.4 and x = 0.5 locate on the fitted curves and the 
crossover point can be seen at x s=s 0.36 which agrees reason- 
ably well with the experimental observation of the crossover 
at x = 0.385 ± 0.016^22 

Since the direct-to-indirect gap crossover is where the con- 
duction band minima at V and X have the same energy 
value, it would be very useful to plot band structures of 
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FIG. 2: (color online) The calculated band gaps of AlzGai-a; As ver- 
sus x by the CPA-MBJ approach. The two solid lines are fitting to 
the data in the two ranges of x. Inset: the calculated DOS for the al- 
loy Alo.36Gao.64As in logarithmic scale as a function of momentum 
k and energy E, revealing a broadened "band structure". 



the Al^Gai-^As alloy. This is however not possible be- 
cause when Al is randomly doped into GaAs, the material 
is disordered and momentum is no longer a good quantum 
number. Nevertheless, one can plot the calculated density 
of states (DOS) as a function of the momentum k and en- 
ergy E as shown in the inset of Fig. 2 which is for the al- 
loy Alo.36Gao.64As. This inset reveals a broadened "band 
structure" of the disordered alloy where the "bands" are no 
longer infinitely sharp lines as that of crystals due to impurity 
broadening. The broadened "bands" trace out the energy min- 
ima and dispersion from which we found that the "conduction 
band" minima at V, X, and L points have essentially the same 
energy value for this alloy, thus confirming that the theoretical 
crossover point is at x ~ 0.36. 

Band offsets of hetero junctions. Having correctly deter- 
mined the band gaps, dispersions and the direct-to-indirect 
gap crossover for the alloy Al^Gai-^As, the band offset of 
the GaAs/AlGaAs heteroj unction can be confidently analysed. 
The valence band offset (VBO) and the conduction band off- 
set (CBO) of a heteroj unction are defined as the difference be- 
tween the energy values of the top of the valence bands and the 
bottom of the conduction bands of the two materials forming 
the junction, respectively. As discussed in the introduction, 
band offsets are of great importance for transport properties 
of the heterojunction. The band offsets can be calculated as: 



VBO, CBO = AE V 



AV 



(2) 



where AE V (AE C ) is defined as the difference between the 
top (bottom) of the valence (conduction) bands of two inde- 
pendent bulk materials that form the heterojunction, AV is the 
lineup of the potential through the heterojunction. When the 
two materials on either side of the heterojunction is extended 
long enough, AV is independent of the orientation of the two 
material interfaces. 
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FIG. 3: (color online) Valence band offset (red dot) and conduc- 
tion band offset (blue square) at different concentration x. The red 
line shows the linear fitting of the VBO. The GaAs/AlzGai-^As 
heterojunctions have the straddling type gap - the valence band 
maximum (VBM) of GaAs is higher, while its conduction band 
minimum (CBM) is lower. Insets: band alignment diagrams of 
GaAs/Alo.sGao.sAs (on the left) and GaAs/AlAs (right). The VBO 
(CBO) values are with respect to the VBM (CBM) of GaAs in units 
ofeV. 



Previously, by supercell method the VBO of the 
GaAs/Al-rGai-zAs heterojunction were calculated at a few 
particular values of x, and several different scalings of VBO 
versus x were reported. For instance Ref. |23] reported VBO 
- 0.564a: - 0.032a: 2 (units eV); Ref. reported VBO ~ 
0.41a; 2 ; and Refill reported VBO ~ 0.17a:. Recently, Ref. IH 
reported a calculation at x = 0.5 by the hybrid density func- 
tional method and obtained VBO=0.26 and CBO=0.42. How- 
ever, first principles analysis of band offsets for the entire 
< x < 1 range is still lacking. 

By the CPA-MBJ method discussed above, we have cal- 
culated VBO and CBO of GaAs/Al^Gai^As heterojunc- 
tion for the entire Al concentration range, results plotted in 
Fig. [3] versus x in intervals of 0.1. The two insets of Figj3] 
plot the band alignment diagrams for x = 0.5 and 1.0. It 
is apparent that VBO is a linear function of x while CBO 
is not. By fitting VBO linearly (solid red line) we obtain 
VBO(x) ~ 0.587a: eV which agrees quite reasonably with 
the experimental scaling 3 of VBO(x) ~ 0.55a: eV. The calcu- 
lated VBO values in Fig. [3] agree very well with the previous 
experimental and theoretical reports, 2 ' 27 namely at x = 0.3, 
the calculated VBO is 0.167 eV, to be compared with the ex- 
perimental value of 0.17 ±0.04 eV.— Concerning CBO, it was 
known experimentally that CBO appears to fit an approximate 
relation CBO~ 0.6AEg where AEg is the band gap differ- 
ence of the two parts of the heterojunction.^^ The linear in- 
crease of the CBO in the direct gap range in Fig. [3] is very 
close to this relation. Near and beyond the direct gap range, 
the calculated CBO decreases with increasing of x, which is 
also in good agreement with the experimental results.— 

Summary. Using the DFT-CPA-MBJ first principles 



approach we have calculated the band offset of the 
GaAs/Al a .Ga 1 _ :I . As heterojunctions for the entire range of the 
Al doping concentration x. The calculated offset scales with 
the Al concentration a; in a linear fashion as VBO{x) ~ 
0.587a;, and the conduction band offset scales with a: in a non- 
linear fashion. In our calculations, the impurity doping is han- 
dled by the CPA while accurate band structures and band gaps 
of the materials are obtained by the MBJ semilocal exchange 
potential. Our atomistic calculations of the band structures of 
GaAs and AlAs crystals, band gaps of the A^Ga^^As al- 
loy, and band offsets of GaAs/Al^Ga^^As heterojunctions 
compare very well with the corresponding experimental re- 



sults. From the calculated density of states, a broadened "band 
structure" at the direct-to-indirect gap crossover point is ob- 
tained at which the alignment of the conduction band minima 
at the T, X and L points occur, also in agreement to the ex- 
perimental observation. 
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